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We propose a robust, unified framework, in wliicli tlie similar baryon and dark matter cosmic 
abundances botli arise from tlie pfiysics of weakly interacting massive particles (WIMPs), with the 
rough quantitative success of the so-called "WIMP miracle". In particular the baryon asymmetry 
arises from the decay of a meta-stable WIMP after its thermal freezeout at or below the weak scale. 
A minimal model and its embedding in R-parity violating (RPV) SUSY are studied as examples. 
The new mechanism saves RPV SUSY from the potential crisis of washing out primordial baryon 
asymmetry. Phenomenological implications for the LHC and precision tests are discussed. 



INTRODUCTION 

The observed dark matter (DM) and baryon abun- 
dances r^DM — 23%, Hb — 4% have long been addressed 
with separate mechanisms at separate scales. The con- 
ventional paradigm for DM theory is the "WIMP mira- 
cle" which gives a striking yet rough guideline for r^DM^ 
thermal relic abundance of a stable WIMP naturally falls 
in the right ballpark of the observed J^dm- The past few 
years have seen rising interest in the intriguing "coinci- 
dence" of JIdm ^ ^b, bringing in the new paradigm of 
"Asymmetric Dark Matter" (ADM)fl. However, ADM' s 
success is at the cost of the WIMP miracle. A unified 
mechanism that can both address the "coincidence" and 
preserve the WIMP miracle would surely be more desir- 
able. Only very recently, a few attempts have been made 
in this direction However, Q is highly sensitive to 
multiple initial conditions; Q proposes a novel baryogen- 
esis triggered by WIMP DM annihilation, but moderate 
adjustment of parameters is required to suppress washout 
effects; Q is also sensitive to washout, and its reliance on 
leptogenesis further restricts working parameters. In this 
paper we explore an alternative baryogenesis mechanism 
with a robust connection to the WIMP miracle and less 
sensitivity to model details. 

Various scenarios addressing the electroweak Hierar- 
chy Problem come with new particles of WIMP type. 
Generically there may be an array of WIMPs, some of 
which are stable, some of which decay promptly, some 
of which have long lifetimes, depending on protection 
from symmetries and mass hierarchies. Although con- 
ventionally, the WIMP miracle only applies to stable 
WIMPs as DM candidates, it has more general appli- 
cation. We consider a meta-stable WIMP that first un- 
dergoes thermal freezeout and later decays in a 0,P^ 
way, triggering baryogenesis The resultant VLb in- 
herits the would-be miracle abundance from the WIMP 
parent up to only moderate suppression from CP asym- 
metry and baryon/ WIMP mass ratio, and thus makes it 
roughly comparable to JIdm of a WIMP DM. A precise 



fit to r^B, i^DM only requires 0(1) adjustment of the dif- 
ferent WIMP parameters, and is insensitive to the precise 
WIMP lifetime. Furthermore when embeded in RPV 
SUSY, this mechanism provides a remedy to a cosmo- 
logical problem there: leading to prompt decays at 
collider typically washes out primordial baryon density 
and calls for baryogenesis below the weak scale. Alter- 
native solutions to this problem 0, 01 are less generic. 
Refs.fs', 'i'l considered low scale baryogenesis in ^SUSY 
to solve the gravitino problem, but the results are sen- 
sitive to details about the infiaton or gravitino. These 
works do not address the WIMP miracle or JIdm — 
"coincidence" . 



GENERAL FORMULATION 



Stage- 1: WIMP freezeout 
A thermal WIMP x freezes out of equilibrium around 
Tf when its thermal annihilation rate Ta — n'^{aAv) 
matches Hubble rate H. This results in the estimate [lo|: 
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which is typically ~ ^"^x- 9 counts the internal degrees 
of freedom of x- 9* counts total degrees of freedom of 
relativistic species. At the end of this stage the co- moving 
density of x is: 
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where s is entropy, g*s = If X is stable, Y^{T{) ~ 

Y^(To), where Tq is today's temperature, and its relic 
density today is: 
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where po = -^o and Sq are the current energy den- 
sity, Hubble rate and entropy, respectively. The second 
line in eq.® manifests the dependence on model pa- 
rameters in the generic case of heavier mediator with 
Now consider two species of WIMPs: 
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Xdm which is stable DM, and xb which decays at time 
T, after freezeout. The observation that eq.dl]) readily 
fits the measured dark matter abundance I^^dm — 23% 
is the well-known "WIMP miracle". In case of xb, 



F^"', acts as the initial condition for later 



baryogenesis, as we now discuss. 
Stage-2: Baryogenesis 

Consider the baryogenesis "parent" xb to have J3f^, 
decay after its freezeout but before BBN, i.e. 1 MeV ^ 
Tbbn < Td < Tf, so that we can treat freezeout and 
baryogenesis as nearly decoupled processes, and retain 
the conventional success of BBN. Solving the Boltzmann 
equationsflO] we get the asymmetric baryon density per 
co-moving volume today YB{Ta ~ 0): 
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where we assume xb decay violates _B by 1 unit, ecp 
is CP asymmetry in xb decay, P-yv is the rate of 
washout processes. Y^"' represents possible pre-existing 
B-asymmetry, which we first assume to be 0. In case of 
weak washout, i.e., Pw < the exponential factor in 
eq.dl]) can be dropped. Then using eqs. p|4|) we obtain: 
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where flZ'^°° is the would-be relic abundance of WIMP 

XB 

Xb in the limit it is stable, given by eq.(|31). ^b given in 
eq.(l5]) is insensitive to the precise lifetime of xs as long 
as it survives thermal freezeout. The observed fJs ~ 4% 
today corresponds to YbIO) = '-^ ^ 10"^°. rJB(O) in 
eq. (HJ takes the form of WIMP miracle, but with an ex- 
tra factor ecp-^^ ^ 10'* — lO^'' for weak scale xb and 
0(1) couplings and phases. Nonetheless as can be seen 
from eq. ^ , the observed ss i can readily arise from 
0(1) difference in masses and couplings associated with 
the two WIMP species xdm and xb- This is our central 
result. 

Note that as long as x decays well before BBN, the 
produced baryons get thermalized efficiently, because 
^px^px ~ r > 7J at Tbbn <^ T < Tew, where 
X can be e^,p,p in the thermal bath. Thus as in 
conventional baryogenesis, the symmetric component of 
baryons is rapidly depleted by thermal annihilation. Di- 
lution/reheating from Xb decay is negligible because at 
Td the energy density of xb is much less than radiation 
density. To see this, recall that today Tq « lO^^eV, 



« lO^. Red-shifting back to Td and using eq.® 
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and sizeable ecp 
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MINIMAL MODEL AND CONSTRAINTS 



We add to the Standard Model (SM) Lagrangian: 
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AC = Xij(j)d,dj + Eixu^ct) + M'^x' + Vii^Ui4> + M|V^ 



ax^S + I3\HYS + M^gS' + h.c. 



(6) 



where H is the SM Higgs, d, u are RH SM quarks, with 
family indices j — 1, 2, 3, is a di-quark scalar with same 
SM gauge charge as u. x, ij} are SM singlet Majorana 
fermions, and 5 is a singlet scalar. x = Xb is the earlier 
WIMP parent for baryogenesis. ^ 1 are our formal 
small parameters leading to long-lived x- They can rep- 
resent a naturally small breaking of a x-parity symmetry 
under which only x is odd. S mediates thermal annihila- 
tion of XX into SM states. The first 3 terms of eq.® give 
rise to collective breaking of U{1)b- Out-of-equilibrium 
decay x — > 0*u is followed by the prompt decay (p ^ dd 
with AB — l,ecp 7^ 0- CP asymmetry ecp in x decay 
comes from the ■(/'-mediated interference between tree- 
level and loop diagrams as shown in Fig.([T|). In the case 
of > 
obtain: 



M^, in close analogy to leptogenesis[llj. we 
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which is non-zero for generic complex couplings. We 
also see that the key to a large ecp is to have yi ~ 0(1) 
for at least one flavor i. Note the analogous ecp from ij; 
decay is 0(£^), with e o y,M^ o in eq.©. 

It is straightforward to incorporate WIMP DM by 
introducing another singlet xdm with analogous interac- 
tions to X) except with £dm = enforced by an exact 
XDM-parity. We will not write out the xdm physics 
explicitly. We next consider various constraints on this 
minimal model. We start with a generic flavor structure, 
and drop family indices in y,e for now. 
Lifetime of x- 

2 

The decay rate of x at T < is Pd — '^g™^ . With 
Tf ~ 100 GeV, our requirement of x decay within 
range Tbbn < Td < Tf leads to the constraint: 
10"" <e< 10~^ 

We next consider potential washout effects. We 
will focus on considering processes involving there 
are analogous diagrams with ^ — ^ Xi bnt they 
give much looser constraints since e <C y ~ 
Early time washout: at T > Aqcd 

As we will see, in this epoch Tw/H decreases with T. 
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FIG. 1. Loop diagrams that interfere with tree-level decay to 
generate ecp 



Thus for each early washout process X, we define 
such that < H ioY T <T^. We require To < 
to have weak washout effect. 
A. Inverse decay udd — >■ via an onshell 0*: 



This gives the constraint: 
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B. AB = 1,2^2 scattering -^u — via ^-exchange: 
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C. Ai3 = 23—^3 scattering udd — >■ udd via on-shell </) 
and V'-exchange. This is effectively 2 — >■ 2 {(t)*u — >■ 
and similarly to case B: 



AS=2.2-i.2 
W 
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D. AS = 2 3 ^ 3,2 ^ 4,4 ^ 2 scattering: udd 
udd via "^-exchange and offshell (p, or ltd — > uddd. 
Taking 3 — >■ 3 for example: 



167r(27r)3 m^m^ 
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(13) 



3^3 /L66g|128^\ m ^^^^ 

where we simplified the expression by taking all masses 

We compare the constraints on Td given in eqs. ([9l 
m] [T2| [T4l) with Tf given in eq. ([T]) where for this model 

p-wave annihilation (ctaw) ^ i^^^^ l^i' "^x < ^i^'l 



0(1) couplings, ~ With non-hierarchical weak 

scale masses of x, V': 0: and 0(1) couplings, we find for 
all washout processes considered, Tw ^ Tf. Therefore 
with Ty) < Tf, early washout is not a concern. Notice 
that potential washout from EW sphaleron is also easily 
avoided since sphaleron shuts off at ^ 100 GeV Tf > 
Td for up to 0(l)TeV. 

After the QCD phase transition, the neutron and pro- 
ton become new effective degrees of freedom to consider. 
n — n oscillation is the typical washout process in this 
era. The general formula for the transition probability is 
0: 
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where Sm is the JS'Majorana mass. The splitting IS.E = 
En — En is in vacuum or in medium where n, n are sym- 
metric, e.g. thermal bath shortly after QCD transition 
when baryons are dominated by the symmetric compo- 
nent. /S.E 3> 5m may occur in an asymmetric medium, 
e.g. the thermal bath close to BBN time or the nucleus 
environment after BBN, which strongly suppresses Pn^n ■ 
In a medium where there is a characteristic time scale r, 
the washout rate can be estimated as 
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Pn^n{T)/T. 



(16) 



Intermediate-time washout: T < ^qcd 
In this epoch n scatters off the thermal background and 
r is set by the mean free path of n, bound by from 
above. In reality both /S.E and r are varying functions in 
this period. To simplify we consider the most "danger- 
ous" limit where /S.E — and r which max- 
imizes washout according to eqs. (ll5ll6p . p^"''"'™ ~ 
{5mfH-\ Requiring r^"''""" < H a± T < Aqcd, 
we find 6m < lO-^SQeV. 
Late-time Washout: T < Tbbn 

After BBN, n is bound in the nucleus. Now the char- 
acteristic time r is set by nuclear time scale which is 
Tnuc (IGeV)-i. In nucleus AE - 100 MeVfl^. Thus 
in this era eq. p5)) becomes approximately: Pn^n ~ ■ 

ft, late Sm^ I p„_ 

► y ^ "(AB)^ / 'nuc • 

late ^ rr c„ J <r i n-22/ 



Thus the washout rate is 



quiring F^^ 



< -ffo, we find 



GeV. 



Current day precision tests: 

n — n oscillation reactor experiments today set a bound 
8m < 6 X 10~'^'^GeV ~ (10^sec)~i which is stronger 
than the washout constraints above. Now we consider 
constraints from 8m on model parameters Ay. In this 
minimal model, Ai^ for (j)didj have to be anti-symmetric 
in z, J . Consequently uddudd operator giving rise to 6m is 
highly suppressed, and Xtj are not effectively constrained 
by rt — 71 oscillation [13|. More relevant constraint comes 
from pp — ^ K'^K'^ decay via higher dimensional ^ op- 
erator, which gives bound A12 ^ 10^^ for m^jm^f, ^ 1 
TeV 13]. As we will show later, when embedding this 
model in natural SUSY where additional fields such as 
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di and related interactions are involved, n — n oscilla- 
tion gives strong bound on A- type couplings. We are 
also constrained by flavor changing neutral currents such 
as Dq — Dq mixing, which gives 2/12/2 ^ 10~^ with TeV 
masses. The large ecp required for baryogenesis may 
bring additional constraints from the neutron EDM. If 
ecp comes from an 0(1) phase in m^, yi, in the minimal 
model where new couplings only involve RH Ui, then the 
contribution involving external quarks vanishes at 1-loop 
for a similar reason as in the SM 14, 151: the dominant 



contribution then arises from the Weinberg operator at 2- 
loop, which still allows phase up to 1 /3 for 0(1) couplings 
and TeV masses An even safer option is to have large 
ecp come from phases Tn^,ei, so that the EDM is safely 
suppressed by ~ jlp- < 10~^^. 

Now we have seen that precision constraints require 
the new couplings to the first two generations of quarks 
need to be suppressed. A simple solution is to consider a 
third-generation dominated pattern where the new fields 
couple mostly to b, t, with CKM-like suppressions to light 
quarks. This choice further strongly suppresses the ear- 
lier washout. 



SUSY INCARNATION AND PHENOMENOLOGY 

Now consider an incarnation of our minimal model 
in "natural SUSY" ^ framework with RPV 
couplings fisj. We promote singlets x ^-nd S to chiral 
superfields which we add to the MSSM. Superpotential 
terms relevant to our setup are: 

W D KjTD.D, + e'xHuHd + VtQHuT + +^i^x^ 
+ fiHuHd + fisS^ + ax^S + f3SHuHd. (17) 

We assume SUSY breaking such that scalar component 
of X and the first two generation squarks are heavy and 
decouple from the low energy spectrum, as in "natural 
SUSY" . The diquark (f) in our minimal model is iden- 
tified with the light tn in superfield T, Majorana 4' is 
identified as a gaugino (Dirac Higgsino mass is not ^). 
In eq. P7|) the terms in the first line ensures and fJf^ 
in X decay, the /z-terms give masses to fermions and also 
induce S — mixing which enables a promising chan- 
nel for LHC search as we will discuss later, and the last 
two trilinear terms involving S provide WIMP annihila- 
tion for X- is a reflection of the e in our non-SUSY 
model, enabling late decay x ^ tt via x ~ Hu mixing. 
Most of our earlier analysis for the non-SUSY model 
directly applies here, except for effects from additional 
fields and interactions. Here gaugino 1/; has both LH 
and RH couplings. Therefore if ecp is from a gaugino, 
the 1-loop neutron EDM with external quarks is non- 
vanishing, but is well suppressed with third generation- 
dominated flavor pattern[l6[. n,n oscillation now con- 
strains A12, A31 < 10"'^, but A23 could be 0(1)|T3], which 



are again naturally satisfied with third generation domi- 
nance. 

RPV ^ natural SUSY is intriguing in both theoret- 
ical and experimental aspects. However, this scenario 
suffers from a cosmological crisis. Assuming an oth- 
erwise successful conventional baryogenesis at or above 
EW scale, RPV strong enough for prompt decays within 
the LHC would typically wash out any primordial B- 
asvmmetrv[20l|. Our SUSY model serves as a robust 
cure to this problem by having baryogenesis below the 
weak scale when all wash-out effects decouple. To see 
the problem clearly, as shown in [18], for a natural stop 
that dominantly decays by ^couplings, A^- > 10"'' is re- 
quired to have prompt decay at collider, i.e. decay length 
L < 1mm. On the other hand, Ay > 10^^ happens to 
be the range where ^scattering such as Hut — > did; can 
efficiently destroy pre-existing B-asymmetry Y^"'* [20| . A 
simple estimate of such wash-out effect can be read off 
by dropping the first term at RHS of eq.([3]). With 



XfjUfT, we find an exponential reduction Yb{T ; 



Y^"'*e 
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LHC Phenomenology 
A promising channel is single resonance production of 
a mostly-singlet heavy scalar admixture of H and S 
which dominantly decays to XX- The production chan- 
nels are the same as for the SM Higgs, except for a mix- 
ing suppression. At 14 TeV LHC run, a Higgs-like boson 
can be produced copiously, even when it is as heavy as 
800 GeV, with say 10% mixing, - 10 fb. The pro- 
duced X must live beyond its freezeout, so its lifetime 

TD >tf (Isec) > 1cm, where T/ < 100 GeV 

so that < O(TeV) is within the LHC reach. Close 
to this bound on tu, x decay leaves a displaced vertex 
inside detector involving t, t. The search can be based on 
dedicated displaced vertex trigger pi, 22|, or triggered on 
two tagging jets in VBF production channel. A challeng- 
ing but exciting further step is to measure {3^ responsi- 
ble for baryogenesis from the charge asymmetry in the tt 
system. 



SUMMARY/OUTLOOK 

We proposed a new mechanism addressing JIdm — 
"coincidence" while preserving the merits of the WIMP 
miracle, presenting a simple example model as well as its 
incarnation in ^natural SUSY. Our basic idea allows for 
further elaborations, e.g. the WIMP parent may decay 
to both asymmetric DM and baryons, or baryogenesis 
may proceed through 3-body deca y, a ccommodating a 
lighter x- In mini-selit SUSY the latter has a 

natural incarnation (26|. On the phenomenology side, our 
mechanism brings the exciting possibility of having the 
cosmological origin of matter being testable at current- 
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day colliders. It is also possible that with improvements 
low energy experiments will be another frontier to test 
the mechanism we proposed. 
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